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ABSTRACT 

We present a hermitian matrix chain representation of the general solution of the Hirota 
bilinear difference equation of three variables. In the large N limit this matrix model 
provides some explicit particular solutions of continuous differential Hirota equation of 
three variables. A relation of this representation to the eigenvalues of transfer matrices 
of 2D quantum integrable models is discussed. 
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1 INTRODUCTION 

The Hirota bilinear equations (HE) [||] provide, may be, the most general view on the 
world of exactly solvable models, from integrable hierarchies of differential and difference 
equations equations, like KdV equations or Toda chains, to, as it was recently shown, 
the transfer matrices of 2D models of statistical mechanics and quantum field theory 

[5], [L|] integrable by the Bethe ansatz (BA) techniques. Indeed, in the last case, it 
was shown that the eigenvalues of transfer matrices must obey the general HE with 3 
discrete variables corresponding to the rapidity, rank and level of a representation with 
rectangular Young tableaux, used in the fusion procedure. The Thermo dynamical Bethe 
Ansatz (TBA) equations discovered by Lee and Yang and widely used in the last years 
for the description of thermal properties and finite size effects of 2D integrable models 
follow almost directly from this HE [Q]. 

HE has its own history in connection with the matrix models (MM). The general 
solitonic solutions of HE [@] bear a striking resemblance with the matrix models with 
logarithmic potentials. The discovery of double scaling limit in the matrix models [0. |9|. 
pTQfl 1 (corresponding to the big size of matrices and a special tuning of potentials) showed 
that the matrix models are closely related to KdV and KP hierarchies of integrable 



differential equations [[ I L]]. The study following it provided even more general examples 



of this correspondence: the multi matrix models before taking any large N limit appeared 



to be related the classical Toda chains [p2|1. Another manifestation of these connections 
are the Schwinger-Dyson equation for MM which can be written in terms of Virasoro 
constraints [0, [14]]] (see [|15[j1 for a modern account of this approach). 



We would also recall one rather mysterious coincidence: in [ ||16|| 1 an open string am- 
plitude for the (1+1)D string with mixed (Dirichlet-Neumann) boundary conditions ap- 
peared to have the same form as the Sine-Gordon S-matrix. 

All this suggests that matrix models could have something to do with the quantum 
2D integrable models and the way from one to another might go through HE. 

In this paper we will propose a matrix chain representation of the general solution 
of Hirota difference equation. We will use for that the analogy between the so called 
Bazhanov-Reshetikhin determinant representation [[17]] of the transfer- matrix eigenval- 



ues of integrable models (which obeys the HE) and the determinant representation of 
matrix chain, in terms of eigenvalues of the matrices. The potentials acting at every site 
of the chain will depend on the eigenvalues of the matrices and on the coordinate of the 
site. Hence the arbitrary potential is a function of 2 variables which is in general enough 
to parameterize any solution of a difference equation of 3 variables. 

The solutions of HE relevant to TBA must obey very special boundary conditions. It 
can include, for example, the condition on the maximal possible size of the fusion repre- 
sentation (reflecting the invariance of the initial model with respect to some continuous 
symmetry), and Lorentz invariance of the spectrum of physical particles emerging on the 
top of physical (dressed) vacuum. 

It is not easy to extract the physical information from HE or the corresponding TBA 
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equations as it is not easy to solve the nonlinear difference or integral equations with 
specific boundary conditions. The results obtained on this way are quite limited: they 
mostly concern the calculations of central charges and dimensions of operators in the 
conformal (ultraviolet) limits and various asymptotic expansions corresponding to high 



energies (see for example [|17|, [18], [HJ]). It even appears to be difficult to reproduce 
the first loop calculation for the asymptotically free models in finite temperature, apart 
from some simplified models [0]), although numerically the TBA equations work quite 
well. Another challenge is to find the planar (large N) limit of such an interesting 2D 
quantum field theory as the principal chiral field. The model was formally integrated by 



the BA approach in [J2T], |2"2fl l. It has been solved rather explicitly in the large N limit in 



case of zero temperature and arbitrarily big external field [ p3|, g4J ], but the attempts to 
generalize it to finite temperatures were not successful. 

On the other hand, the general solution of HE follows from its integrability and can 
be represented in terms of a r-function. In case of a general difference HE it coincides 
with the determinant representation of Bazhanov-Reshetikhin. It might be a good idea 
to use this representation, and hence our matrix chain representation, to get some hand 
on HE and TBA. Of course, for a finite rank N of matrices (and hence finite dimensional 
groups of symmetry of corresponding models) our representation hardly could offer some 
breakthrough. But for big N we can try to use the machinery of the matrix models (like 
orthogonal polynomials, character expansions and various saddle point techniques) to 
calculate the corresponding infinite determinant. 

We were not able to find any solutions of HE satisfying correct TBA-like boundary 
conditions. A natural way to impose these boundary conditions is the most important 
drawback of our representation. Leaving it to future studies we propose here some par- 
ticular solutions of continuous (differential) Hirota equations which describe our matrix 
chain in the large N limit. They correspond to a particular choice of parameters (poten- 
tials) of the chain. 

In the next section we will briefly review how the Hirota equation is connected with 
the integrable models of the 2D quantum field theory. The continuous differential version 
of them will be presented. 

In section 3 we will propose the matrix chain representation of the difference Hirota 
equation based on the Bazhanov-Reshetikhin determinant representation of the fusion 
rules for transfer-matrices. 

In section 4 we will present some particular examples of solution of the continuous 
differential HE, given by the one matrix model and the matrix oscillator with the specific 
boundary conditions. 

In section 5 we will sketch out a general solution of the differential HE for an arbitrary 
time dependent matrix chain potential and consider a more explicit solution for the 
particular case of time independent potential. 

The last section will be devoted to conclusions and prospects. 
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2 TBA, FUSION RULES 

AND HIROTA DIFFERENCE EQUATION 

To set a more physical background for our construction let us briefly review how the HE 
appears from the Bethe ansatz. We will mostly follow in this section the framework and 
the notations of [0,0] . 

The transfer-matrix of an integrable 2D model with periodic boundary conditions 
depends (apart from a number of fixed parameters, like volume, temperature or the 
anisotropy q) on three variables: rapidity u = id, rank ("color") a and level ("string" 
length) s. The variables a and s have the meaning of a representation of elementary spins 
filling the bare vacuum of the model, given by the rectangular Young tableau of the size 
a x s. The corresponding transfer matrix is called Tg(u). 

The integrability imposes the commutativity of transfer-matrices for different values 
of all three variables playing thus the role of spectral parameters: 

[t>),t;V)] = o (i) 

It follows from ([!]) that we can always work with the eigenvalues T s a (w) instead of the 
transfer-matrix itself and view them as usual functions. 

The transfer-matrices, as well as their eigenvalues, obey a set of relations known 
as fusion rules, originally found as the relations between S-matrices for particles with 
different spins in integrable QFT. They can be summarized in the so called Bazhanov- 
Reshetikhin formula [ ||17|| 1 (BR) presenting the function of three variables T"(w) in terms 
of the function of only two variables T}(u): 

T s » = dehKijKaT^iu + i + j + a) (2) 

Actually, there exists a more general BR formula, expressing the transfer-matrix 
eigenvalue of any skew representation h/h' through T}(u): 

Th/h>{u) = det 1 < i j< a Tl._ h ,{u + hi + h'j) (3) 

where hi = rrii+a—i and h\ = m^+a—i are the so called shifted highest weight components 
of two representations R and R' of GL(N) characterized by the usual highest weight 
components R = (m\,...m a ) and R' = (m' 1; ...m' a ), so that they obey the inequalities 
hi < hi-i, h\ < h'i and also h\ < hi. But the transfer-matrices with rectangular Young 
tableaux play an exceptional role since they obey a closed set of fusion rules given by the 
difference Hirota equation: 

T> + 1)T S > - 1) - T^T^u) = T«+\u)Tr\u) (4) 

It follows directly from (0) in virtue of the Jacobi identity for determinants. It contains 
little information since it is true for any function of two variables T]{u) in (fj). To 
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specify it further to some particular integrable model we have to impose some boundary 
conditions on solutions of the eq. @. 

One of these conditions specifies the group of symmetry of the model. To make it, 
say, SU(N) (or A^^i in terms of underlying algebra) we put: 

T?(u) = 0, for a < and a> N (5) 

It is not enough since it leaves us with an infinite discrete set of possible solutions 
(like in quantum mechanics, fixing the boundary conditions on a wave function we are 
still left with infinitely many wave functions corresponding to different energy levels). 
We have to specify some analytical properties of solutions. 

There are two ways to do it in the case of BA. 

One is related to the so called bare BA where one specifies T s °(u) and T^(u) to be 
some given polynomials in the variable u whose zeroes specify completely a model, where 
as the functions T®(u), for 1 < a < N — 1 are polynomials whose zeros we have to 
find. The details of the analyticity conditions for the bare BA can be found for example 
in (I). 

Another way to fix analytical properties corresponds to the dressed BA where the 
elementary excitations are already the real physical particles. To precise them let us 
derive from @ the TBA equations. For that we introduce the function: 

ya, u) = TUU)TUU) 

which, in virtue of (§), satisfies the equation sometimes called Y-system [||]: 

Y?(u + l)Y?(u - 1) [1 + + Y a U(u)] 

Y s a + l (u)Y s a -\u) ~ [l + Y s a + 1 (u)}[l + y/" 1 ^)] 

Note that this system is symmetric under the change: Y — > Y~ x , a — > s, s — > a 
(rank- level duality). 

To make it a little bit more symmetric let us introduce the functions: 

c/» = i + y s », [/» = !+ * (8) 

Y s a {u) 



(7) 



Then the eq. (0) can be rewritten as 

U?(u + l)U^(u - 1) U*(u + l)U?(u - 1) 



UUi (u) UUi (u) (it) u-- 1 («) 

Taking logarithm of both sides of this equation and applying the operator 



(9) 



oo ] 

-co ^shSf* -0'f" (10) 
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where 6 = iu we obtain: 



C s ,{0) * In (1 + Y s %e)) = C aa \6) * In (1 + — (11) 
where we introduced the so called "baxterized" Cartan matrices: 

and similarly for C ss >. The * sign defines the usual convolution operation: f(8) * g(6) = 
I^de'f{6-e')g{9'). 

Let us now act on both sides of eq. ([□]) by the operator inverse to (|12"D: 

Aaa' 2 f°° , f n\ ,smh\(N - max(aa'))p] sm\i\min(aa')p] . . 
a n = Cj(9) = -J Q dpoos(pe)coth(p) A_^L ^_Jfl (13) 

Note that this inverse is respecting the boundary conditions restricting the values of a, a' 
to 1 < a, a' < N - 1. 

Note also that the operator C aa has zero modes: 

C aa \6) * a{s)m sm{a'irk/N) cosh(#vr£;/iV) = (14) 

for any integer k and any function cr(s). So in acting by ( fl3|) on both sides of (|TT1) we 
might be obliged to add one of zero modes. The choice of zero mode and the function a(s) 
defines completely the boundary conditions and hence the model. If we want to respect 
the 1+1 dimensional relativistic invariance we can add the zero mode with k — 1, since 
only it will lead to the relativistic spectrum of energies elementary excitations (which are 
described by this zero mode) of a type: 

a(s)m sin(a7r/iV) cosh(6>7r/iV) (15) 

It gives a typical mass and energy spectrum of physical particles for integrable relativistic 
models of 2D QFT. The choice of a(s) and the range of s define a particular relativistic 
model. For example, for a(s) = 6 8) i, s > 1 corresponds to the chiral Gross- Neveu model, 
whereas a(s) = <5 Sj o, — oo < s < oo corresponds to the principal chiral field (PCF) with 
the SU(N) symmetry. 

With all these settings the final TBA equation (or similar equations for the ground 
state of the finite length system with the periodic boundary conditions) takes a familiar 
form (say, for the PCF): 

A aa '(6) * C ss ,{e) * In (1 + Y£( u )) - In (1 + — = 5 sfi m sm(arc/N) cosh(^vr/iV) (16) 

where e(x,y,r) = logY^w) plays the role of the energy density of the excitations char- 
acterized by the rank a and level s. 
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At the end of this section let us comment on the large N limit of the TBA equations. 
It is not obvious to us how to simplify the eq. (0) in this limit, but the difference 
equation equation (0) after introducing the rescaled continuous variables 

r = u/N, rj = s/N, u = a/N (17) 

becomes a second order differential equation for the quantity e(r), u, r) = logF s a (u) Q 

(Ql _ d l) e = (Ql _ ^) i n (i + exp e ) (18) 

This is an integrable classical equation, as it is a consequence of the general HE and 
the determinant representation (fj) for its solution, although the determinant becomes 
functional in the large N limit. 

Let us give also another form of this equation in terms of new variables l,m and 
a ( and the corresponding rescaled variables) A, fi and v defined through the original 
variables u, s and a as follows: 

. I u + s + a m u — s + a a 

A = — = , u = — = , v = — 19 

N 2N N 2N N K ' 

Note that A and /i play the role of "light-cone" variables. 

In their terms we can represent the function e = In Y in the large N limit as 

e = -ln =4 (20) 

where ^ 

F = h'mjv^oo— logT. (21) 



The continuous HE flT8| ) for this function looks in the new variables as 

d x d, ln(exp[-i^J - 1) = d v {d v -dx- d,)Fl\ (22) 
This form of HE will be useful for our matrix model representation of its solution. 

3 MATRIX MODEL CHAIN AS SOLUTION OF 
THE BILINEAR DIFFERENCE 
HIROTA EQUATION 

In this section we propose to parametrize the general solution of HE (|j) by means of 
so called matrix chain integral - a matrix model widely used and investigated in the 
literature (see [^5f] for the details). 



! I thank P. Zinn- Justin for this comment 
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Let us define the following Green's function: 
K a (T, 0(0), 0(T)) = / n d a2 m exp tr[J2 4>(k)(f)(k + 1) - £ V(k, fa)] (23) 

fc=l fc=0 fc=0 

where <p{k) - are a x a hermitian matrices, each corresponding to its site k of the chain, 
with matrix elements 0j 7 (fc), i,j = 1, 2, • • ■ , N, k = 0, 1, 2, • • ■ , T. Two matrices at the 
ends of the chain are fixed. 

The continuous analogue of the chain would be just the matrix quantum mechanical 
Green's function on the interval of time (0,T): 

K a (T, 0(0), $(T)) = J D a2 4>(t) exp -tr J* dt[^ 2 + V(t, 0)] (24) 

with the end point values of the matrix fields also fixed. The matrix quantum mechanics 
was introduced and solved in (p6|). 
Let us now define a new quantity: 

Z a (l,m) = [ d a2 m [ rf a2 0(T)(det0(O))X(T,0(O),0(T))(det0(T)) m (25) 



where K could be any of both Green's functions (^3|) or (j24|). We have added two 
logarithmic potentials at the ends of the chain to describe the dependence on / and m 
introduced by (|19|). 

We claim that the function of three discrete variables 

obeys the HE ([|). More than that, it gives the most general solution of HE parameterized 
by the function of 2 variables - the potential V(t,x). 

The proof goes as follows. If we start, say, from the discrete version we can diagonalize 
each of the matrices in the chain by the unitary rotation: 



J2(^MZa)k(tt a ) kj (27) 



and integrate over the relative "angles" between two consecutive matrices in the chain 
by means of the Itzykson-Zuber-Harish-Chandra formula. This concerns only the first 
term in the exponent in the r.h.s of (|23|). The potentials, including two determinants at 
the ends, depend only on the eigenvalues. The overall result after the angular integration 
will be (see (|23|) for the details of this calculation): 

ZJl,m)= det / dp I ' dqK(T,p,q)p l+l - l q m+] - 1 (28) 

l<i,j<aJ J 

where 



K(T, p,q)= [ if dz{k) exp{ £ z(k)z(k + 1) - £ V(k, z{k))} (29) 

k=l k=0 k=0 
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for the discrete chain, or 



K(T,p,q) = J Dz(t)ex P -{J Q dt[j + V(t,z)]} (30) 

for the continuous quantum mechanics, where z(0) = p, z(T) = q. 

Now we see that due to the determinant representation Q2"B|) the function T°(u) defined 
by the eq. (|26|) has the same determinant form as the BR formula @ and hence it obeys 
the HE (g), if we identify 

Tliu) = J dp J dqK(T,p,q)p^q^ (31) 

The formula (^Bj) generally defines an arbitrary function of two variables s and u. It is 
clear from the fact that it is just the Mellin transform of an arbitrary function ( p0[ ) (or 
(|29"D) in two variables, which is in our case the Greens function of a quantum mechanical 
particle in an arbitrary time and space dependent potential. This potential obviously 
gives enough of freedom to define K(T,p,q) as an arbitrary function of two variables p 
and q. So we proved (or at least made rather obvious) our statement about the generality 
of this representation of solution of the HE. 

We can provide a more general matrix integral giving the parameterization of the 
most general BR formula @: 



T h/h . = J d a 0o y d a <hX[h]{<h)K*{TMi<h)Xw(<h) (32) 

where X[/i](0) is the GL(a) character of the representation characterized by the highest 
weight [h\. It can be easily proved when written in terms of eigenvalues with the use of 
the Weyl formula for characters. 

For any finite N all this seems to be on the edge of triviality: we just defined in a 
sophisticated way an arbitrary function of two variables and built from it the necessary 
determinant. Naturally, we don't expect this representation to be of big use for a finite 
N. The boundary conditions of TBA will be as difficult to satisfy as before. What is 
our major hope is the large N limit of this matrix model which should correspond to the 
large N limit of the integrable models of the type of principal chiral field. In this case the 
determinant in the BR formula is essentially functional, and the matrix models give a 
rich variety of methods for the calculation of such determinants. That why our strategy 
will be the following: we investigate the matrix integrals of the type fl25|) in the large 
N limit for various potentials and look for physically interesting regimes. Things might 
become much more universal in the large N limit, and it could exist a classification of 
interesting regimes, like it was done for the multi-critical points in the matrix models. 
This paper represents of course only a few modest steps in this direction. 

Let us make an important remark concerning the large N limit of the representation 
): the Y variable introduced in the previous section and presented by the formula (EDI) 



S 



in the large N limit (with the matrix chain partition function Z instead of the transfer- 
matrix T) obeys the differential equation ( |l8f) or, in new variables fll~9|) , (0). It is almost 
clear from our definitions and it will be demonstrated in the following sections. 

Another more formal but interesting application of this method could be the search 
for new solutions in the integrable equations of the type fllBp . In the next section we will 
show some particular examples of it. 

4 EXAMPLES OF SOLUTIONS 

OF THE CONTINUOUS HIROTA EQUATION 

We will demonstrate here on some examples limited to particular choices of the matrix 
chain potentials how this relation between the HE and MM works. 

4.1 One matrix model and GL(n) character 

Let us start from the one matrix model partition function with an extra logarithmic 
potential 

Z a (l) = f d a2 (f)(det(f)yexp-NtrV(<j)) (33) 



which, after going to the eigenvalue representation, becomes 

Z a (l)= det [ dpp l+i+j - 2 exp-NV(x) (34) 

l<i ,j<a J 



l<i ,j<a 

It is well know that if one chooses: 

a 

exp-V(x) = H(b k -x)- 1 (35) 

k=l 

and performs the integral in (|34]) along the contour encircling all these poles one will 
identify this partition function with the GL(a) character x?{b) of the a x / rectangular 
Young tableau given by the Weyl determinant formula. So, much of our next formulas 
is valid for the characters as well. 

It is easy to see from the Jacobi identity for determinants that the function 

t a (l) = Z a (a - I) (36) 

satisfies a simplified version of the general HE ([|): 

t a {l + l)t a (l - 1) - t a+1 {l)t a -\l) = [t a {l)] 2 (37) 

Introducing the variable 
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we obtain from (R7I) a simplified version of the general Y-system (|7]): 

Y a {i + \)Y a {i - 1) [i + Y a+l {i)}\\ + r-^)] 



(39) 



If we go to the large N limit and introduce the continuous variables (|19|) we obtain from 



the previous equation a simplified version of the differential HE ([181 ) : 

«9>gF = «9>g(l + r) (40) 

which is again integrable, as it is obvious from the above determinant formulas. In the 
last equations we changed a bit the definition of the function Y rescaling the variables 
byl/iV. 

Take first a simple example of the potential V(x) = \x\. The direct calculation of 
gives: 

4(0 = f[(l-a + k)\(k-l)\ (41) 
fc=i 



which yields Y a {l) as 
or, in the large N limit: 



YJX) = ^ (43) 



which perfectly satisfies the eq. ([40]). 

To find a general (up to some comments which will follow) solution of the eq. ( f4"0"|) 
we just have to apply the well known formulas for the saddle point approximation in the 
one matrix model with the potential which is now V(x) + (A — v) logx. Omitting the 
standard calculations (see for example [0]) we give the result: the function Y U (X) obeys 
the following system of equations on Y and an intermediate variable S: 

•i du V'(VYu + 3) X-v 



I 



1 7T y/T^v? VS 2 - 4Y 

1 du (Wu + S)V'(VYu + S) 



(44) 



i vr y/\ 



u 



2 



X + u (45) 



So we have obtained the solution of the eq. (|40|) in terms of a system of ordinary 
equations. For example, for a polynomial potential the equations will become algebraic. 
In general they are functional. 

Not every potential is compatible with this solution. We restricted ourselves to the 
so called one cut solution implying the existence of one classically stable well in the 
potential. This restricts our solution to some parametrically general but still limited class 
of solutions of (0). The generalization to the multi-cut solution which is straightforward 
should in principle yield the most general solution of ([101). 

The equations ( P| - ^5D look like the characteristics method of solution of the eq. (|40|) . 
In the next sections we shall see to what extent we can generalize it to the full differential 

he ran. 
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4.2 Gaussian Chain Solution of HE 



Now we shell consider another particular example of solution of the HE (f|) by restrict- 
ing all the potentials in the matrix chain to be Gaussian. Then all the integrals over 
0(1), ...,0(T) in ( p3|) can be easily performed and we are left in (p5|) with the following 
two matrix integral over the endpoint variables: 

Z a (l, m)= J d a2 <f) J cf V (det <P) 1 exp -tr(^0 2 + ^0' 2 - c#')(det 0') m (46) 



If we do the same with the continuous quantum mechanical integral (|24|) we arrive 
(up to a trivial coefficient typical for the Green's function of the harmonic oscillator) at 
the same two matrix integral, with c = 1/ cosh(c<jT), where u is the frequency of the 
corresponding oscillator. 

An important comment is in order. Although this partition function satisfies the HE 
for three variables it should a little bit modified for finite N: note that Z a (l,m) in (f46f ) 
after passing to eigenvalues splits into the product of two determinants corresponding to 
i, j both even or both odd (the matrix elements with different parities of i, j are zero, see 
the formula (p8|) with the gaussian kernel). Each of these determinants satisfies the same 
difference HE with the shifts of discrete variables by ±2 and not by ±1. The continuous 
(large N) version of HE will be the same as before. 

This two-matrix model has the only complication with respect to the ordinary one, 
containing usually only polynomial potentials: its potentials contain logarithmic parts, 
like in the well know one matrix Penner model. To solve it the ordinary method of 
orthogonal polynomials does not look convenient. We propose here another, rather pow- 
erful method worked out in a series of papers [ p6| , [2T| , [28| , [29| , |30| , |3Tf l and capable to 



solve some even more sophisticated models than the present one. 

First we perform the integral over the relative "angles" of two matrices in ([|6|) by 
means of the character expansion [ [^2| j: 

/ (dn) u(N) ex P [c tr(<pn + <p'n)} = £ c E fc ^ °(° 1)72 (47) 



0<h a <--<h 1 <oo 



We dropped here some unessential overall coefficient. 

Plugging this formula into the eq. ([R]) we encounter two identical independent Gaus- 
sian integrals over and 0' with the characters as pre-exponentials. These integrals can 
be calculated (they slightly generalize the similar integrals appearing in [R53L 27, [28|, |29| 
to the case of l,m^ 0). The result (again up to some unessential factor) is: 

/cfV-^Mtdet,)' = IMfclV w (48) 

where we denote by h e {h°) the even(odd) highest weights whose numbers should be equal. 
A(/i) is the Van-der-Monde determinant of Ks. We chose here / to be even; for / odd one 
only has to exchange h e and h° in (f48|). 
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Putting all this together we obtain for fllED a representation in terms of the multiple 
sum over h's (we dropped a h- independent coefficient): 



Z a (l,m) 



£ A 2 (h e )A 2 (h )<£*M +h $ x 



(49) 



{h e ,h°} 



X- 



n<(/if + /-!)!! n,(/i° + 0" Mhf + m - 1)!! IL(^ + m)!! 



26 



llrt 

We chose here l,m to be both even. Different parity for them is forbidden by the sym- 
metry — > —0 of (^). This will not be important in the large N limit. As we see, the 
sums over h e and h° are decoupled and can be calculated independently. 

The method of calculation of these multiple sums in the large N limit was proposed 
and further elaborated in [ |43| , |27| , [28| , |29| 1 and is based on the saddle point 



in 

approximation of this sum. One introduces the resolvent function of shifted highest 
weights: 

H{h) = tA (50) 

In what follows we change h by h/N (since the highest weights are supposed to be of 
the order N in the large N limit). So, in the large N limit: 



H(h) 



'to™ 



h-hf 



H+(h) ±inp(h) 



(51) 



where H + (h) is the symmetric part of the function H(h) on the cut defined by the 
distribution of h's and p(h) is the density of h's along this cut. In the large N limit we 
can calculate the multiple sum by the saddle point method. The saddle point condition 
defines the most probable Young tableau shaped by the density p(h): 



4 d dh> 

% h-h> 



In 



'c 2 (h + X)(h + fj,)' 



(52) 



2 ~ \ h 2 J 

One has to remember that a part of the most probable Young tableau is in general 
empty (some of the highest weight components m& are equal to zero, see [[2^, 27\] for the 
details). So, the function p{h) is equal to one on the interval (0, b) and to some nontrivial 
function p(h) on the interval (b,d). This yields, instead of fl5"2"D, the equation: 



r dh! ^ 
t m h~h< 



In 



'c 2 (h + X)(h + p)' 



(h - by 



(53) 



This linear integral equation has a one-cut solution: 



H(h) =ln c'\d-b)h (54) 



2 
1 



In 



(b + d + 2X)h - (b + d)X - 2db + 2J(d + X)(b + X)(h -d)(h- b) 



In 



[b + d + 2p)h -{b + d)p - 2db + 2yJ(d + p)(b + p)(h - d)(h - b) 
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To fix d and b we should recall the asymptotic of H(h) with respect to large h: 

H{h) = u/h+ < h > /h 2 + 0(l//i 3 ) (55) 

following from ([)!]). Here < h > is the average shifted highest weight in the most probable 
Young tableau. Expanding H(h) in 1/h up to the terms 0(1/ h 2 ) we obtain a system of 
equations defining d and b: 



2X + d + b + 2J(d + A) (6 + A) 



2n + d + b + 2J(d + n)(b + (i) 



c- 2 (d-b) 2 (56) 



'(d + A) (6 + A) + yj(d + fi)(b + n) = 2v + A + // (57) 
From (|50| ) we deduce the following formula for the solution of the continuous HE (22): 



F'L = 7 « ln(/i + A), ln(/i + /i) » 



4 



(5? 



where by << A, B » we denoted the connected average of any two /i-dependent func- 
tions A and B. Note that this average has a finite large N limit, as it should be. 

This solution can be brought into a more explicit form: the explicite formula for 



such correlators in the one matrix models was given in [[34, [35]] for W(z,z) =<< 
i i 

z-h)l (z'-h) 



>>: 



W(z,z') 



(z-d){z'-b) 



+ 



[z'-d)(z-b) 



\ (z' -d)(z-b) \(z- d){z> -b) (z- z'Y 



(59) 



Integrating it in z and z' and putting z = —A, z' = —\i we obtain a rather explicit 
solution of differential HE (E2|) 



1 



K» = o [9(\ b, d)g{fi, d, b) + g(X, d, b)g{n, b, d) + ln(A + //)] 



(60) 



where 



g(z, 6, d) = ^J( z + d){z + b) + 



d — b z + d 



z + b 



ln[2J(z + d)(z + b) -2z-d-b] (61) 



and d,b are defined by the eqs. (W^i)- 

In the next section we will give the solution of (|T8| ) in the case of a general time- 
independent potential V(x). We will also reduce the search for the most general solution 
of continuous Hirota equation (|T^) defined by the time dependent potential V(x,t) to a 
simpler problem. 
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5 COLLECTIVE FIELD METHOD 
FOR THE DIFFERENTIAL HE 



Now let us briefly describe how our method works in the case of a the general potential in 
the matrix quantum mechanics defined by fl2~4|) . In this case we can apply the collective 
coordinate method of Jevicki and Sakita [JPJ] which is valid in the large N limit and 
can be applied to the non-stationary saddle point solutions which are needed in our case. 
The details of this approach can be found in [[44, ^2], [45], [37]]]. We will use the results of 
it and apply them to our case. 

In terms of this method the effective action for (PI ) can be written for the density 
p(x, t) of eigenvalues Xk and its conjugate momentum P(x, t) developing in time as: 

S eff [p, P} = Jdx jf dt[pP - -pP' 2 + yp 3 + pV{x, t)\ (62) 

reflecting our specific boundary conditions at the ends of the interval (0, T) following 
from ([25|). 

The equations of motion corresponding to the action (|5"2"D are: 



p + d x (P'p) = (63) 

P+\P' 2 = ^P 2 + V(x } t) (64) 

Differentiating the second one in x and defining the function f(x, t) = P' + inp we rewrite 
this system as only one forced Hopf equation on this complex function: 

d t f + fd x f = -V(x,t) (65) 

We have to impose at any moment, say, at t = the normalization condition 

dxp{x, 0) = v (66) 



Then it will be true at any t due to the condition 

We have excluded the end points of the interval (0, T) in the last equation. The 
logarithmic potentials and the two Van-der-Monde determinants left at the ends of the 
interval can be taken into account as the boundary conditions: 

Re fix, 0) = - + 4 dx' ^H^Z, 6(0) < x < d(0) (67) 

X Jb(0) x — X 

Ref(x, T) = - £^ dx' b(T) < x < d(T) (68) 

or, introducing the resolvent: R±(x,t) = j£$ dx' p ^_'*} ^ ^ p(x , t) , 

f(x,Q) = - + R+(x,0) (69) 
x 
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f(x,T) = -^-R4x,T) (70) 

x 

Hence we reduced the problem of the virtually general (since it is parameterized by 
a general potential depending on 2 variables) solution of the continuous HE (f|) on the 
function of 3 variables u, A et /i to the solution of the differential equation of the first order 
(|65"D (assuming that the complex function / is analytical in their variables) supplemented 
by the boundary conditions (|67| - |68D . The variables v, A et // appear here only as fixed 



parameters. 

It is still quite complicated, although simpler than the original problem and, may 
be, physically more transparent. We don't know how we could simplify it further in 
general case. So let us consider an interesting particular example of the time independent 
potential V(x). In that case the forced Hopf equation 

dtf + fd x f = -V\x) (71) 

becomes completely integrable by the characteristics method 0. 

We will choose A = p, which does not look as restriction if we assume the analyticity 
in A and /i Then due to the time reversal symmetry we have f(x,T) = —f(x, 0). So we 
can say that also p(x,T) = p(x,0). Hence only one of these two boundary conditions is 
independent. 

The first part of the problem is to find the solution of eq. ( |7T|) with fixed endpoint 
density p(x, T) = p(x, 0). The result can be formulated as the following equation on the 
functions already at the endpoints 

T = I 01 '" , dV (72) 
J' ^(soto - V(y)) 

where G(x) = xt{9t( x )) an d xx(g) is the function to be found by solution of 

g(x T ) = ^f 2 (x T ,T) + V(x T ) (73) 
Note that the function G(x) = xo(go(x)) defined by the solution of 

9 = lf 2 (x ,0) + V(x ) (74) 
is the functional inverse of the function G(x) itself which gives the equation of A. Matytsin 

my- 

G(G{x)) = x (75) 

Although the dynamics of the forced Hopf equation is summarized by the relation (|72|) 
the last equation leads to a strong constraint on the analytical structure of the function 
G(x). 

4 I am grateful to A. Matytsin for the explanation of this method and its application to the forced 
Hopf equation 
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Once we found the functional f p ( x ,o)(%, 0) as the solution of eqs. 
much it with our boundary conditions: 

A r d (°) o(x' 0) 
Ref p{X:0) (x,0) = - + 4 dx'^L^L^ b(0)<x<d(0) (76) 
X Jb(0) x — x' 

With a given V(x) this defines the end-point density p{x, 0) which is the only non- 
trivial information we need to find the quantity (^). Probably it is convenient to repre- 
sent this boundary condition as a condition on the large x asymptotic of f(x, 0): 

/(z,0)-w,^^ (77) 

We thus reduced the solution of the continuous HE ( p2|) to some simpler functional 
problem in a particular but rather representative case of the time independent potential. 
The solution with the harmonic oscillator potential obtained in the previous section 
should be also reproducible by this method. 

We can also use these equation to produce more explicit solutions of the HE ( p2|) by 
the method proposed in [[[37]]: one chooses two conjugated roots G(x) and G{x) of an 
algebraic equation x{G) = 0, where x(G) is some polynomial. These two roots satisfy by 
construction the eq. Q75p. Then one has to plug them into ( [72p and solve it as an integral 
equation for V(x). Of course the choice should be limited by the boundary conditions 
©or®. 

It would be interesting to analyze the case of the inverted oscillator potential corre- 
sponding to the 1+1 dimensional non-critical string theory. But this question lies beyond 
the scope of this paper. 



72 -uE) we have to 



6 CONCLUSIONS AND PROSPECTS 

In this paper we proposed a matrix model representation for the solution of the general 
difference Hirota equation. For its continuous analog of the differential HE of second order 
on three variables the solution is represented by the large N limit of the corresponding 
matrix chain. It gives an effective framework for solving the continuous HE in a rather 
explicit way, at least for some particular cases. In particular, the problem can be reduced 
to the forced Hopf equation with specific boundary conditions. 

Many things remain to be understood. First of all it is not clear how to find the 
solutions satisfying the boundary conditions of various 1+1 dimensional quantum field 
theories solvable by Bethe ansatz. Especially how to choose the matrix potentials to get 
the relativistic spectrum for the physical particles and to fix the symmetries of original 
models. Another question: is there some physical interpretation in terms of these in- 
tegrable theories of the time variable t and of the eigenvalue variable x in our matrix 
representation similar to the non-critical (1+1) dimensional string theory (where these 
variables describe the target space of the string [ P2| 0])? It might be for example that 
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the time can be considered as the physical space dimension of the corresponding inte- 
grable theory. But it remains to be proved. We hope that the formalism proposed here 
at least sets a convenient framework for the attempts to find the physically interesting 
solutions. 

A more formal use of our method might be the search for new solitonic solutions for 
the well known integrable equations. For example, the Toda equation d^F = d\d^ exp F 
is just a particular limit of the continuous HE (|i8l) . To our knowledge, the solitonic 
solutions to this equation are not yet found. 

Another interesting question is how the double scaling limit in matrix models is related 
to HE? For example, how to find the corresponding solution for the inverted harmonic 
oscillator giving the description of the (1+1) dimensional string field theory. To answer 
this question as well as many others we have to learn how to deal with the non-stationary 
forced Hopf equation with our specific boundary conditions. The methods worked out in 
the papers [ |3"7| , |3"8| , |3l| f4"0"| ] could be useful for that. 
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